Bulk photons in Asymmetrically Warped Spacetimes 
and Non-trivial Vacuum Refractive Index 

. . . K. Farakos^ * N. E. Mavromatos'* ^and P. Pasipoularides^ ^ 

O ■ "^King's College London, University of London 

O' 



o 

(N 



X 



Department of Physics, Strand WC2R 2LS, London, U.K. 



^ '. Department of Physics, National Technical University of Athens 



Zografou Campus, 157 80 Athens, Greece 



Abstract 



rr^. We consider asymmetrically warped brane models, or equivalently brane models 

r~| I where the background metric is characterized by different time and space warp factors. 

The main feature of these models is that 4D Lorentz symmetry is violated for fields 
fS| ■ which propagate in the bulk, such as gravitons. In this paper we examine the case of 

^^ I bulk photons in asymmetrically warped brane models. Although our results are general, 

r^^ I we examine here two specific but characteristic solutions: 1) AdS-Schwarzschild 5D 

OO ■ Black Hole solution and 2) AdS-Reissner Nordstrom 5D Black Hole solution. We show 

^— ^ , that the standard Lorentz invariant dispersion relation for 4D photons is corrected 

C^^ I by nonlinear terms which lead to an Energy-dependent speed of light. Specifically, 

^"^ ■ we obtain a sub-luminous Energy-dependent refractive index of the form ne//(u;) = 

O ■ 1 + cg 1^^, where co is the energy of the photon, and the factor cg is always positive 

and depends on the free parameters of the model. Finally, comparing the results with 
recent data from the MAGIC Telescope, claiming a delayed arrival of photons from the 
Active Galactic Nucleus of MkSOl, we impose concrete restrictions to the two sets of 
5^ \ models examined in this work. We shall also discuss briefly other possible astrophysical 

constraints on our models. 
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1 Introduction and Motivation 

Theorists, in an attempt to solve the hierarchy problem, invented new string theory models 
with relatively large extra dimensions. The early realization [1] that the string scale is an 
arbitrary parameter, and can be as low as a TeV scale, lead naturally to the consideration of 
models with large extra dimensions [2l [3], introducing the so-called brane-world models, in 
which the standard model particles are assumed to be localized in a 3D brane (our world), 
embedded in a multi-dimensional manifold (bulk). Subsequently, models in which the bulk 
space time is warped have been proposed IH [5j . In such models, the extra dimensions could 
be: (i) finite, if a second parallel brane world lies at a finite bulk distance from our world ^, 
thus providing a new hierarchy of mass scales, or: (ii) infinite, if our world is viewed as an 
isolated brane, embedded in an (infinite) bulk space [Sj. In fact, it is the presence of such 
warp factors that provides [Ij in the case (i) a "resolution" to the hierarchy problem. The 
important problem of stabilization of the extra dimension, in the case of two brane models, 
has been considered in Refs. [6]. 

Such brane or string models lead to no apparent contradictions with the present-day 
observations regarding four dimensional physics [7j, although many falsifiable predictions 
can be made of relevance to either particle physics colliders, such as LHC and future linear 
colliders, or astrophysical/cosmological measurements. In the standard brane world scenario, 
only gravitons are allowed to propagate in the bulk, contrary to the standard model particles, 
which are constrained to lie on the brane. We call this the Standard String /Brane- World 
scenario (SSBW). 

Beyond the SSBW scenario, there are models in which all or some of the standard model 
particles, can live in the bulk. The original scenario of universal extra dimensions, where 
standard model particles possess KK excitations, was presented for compactified dimensions. 



R < 10"^'^ cm, in [8]; the connection with brane worlds of Randall- Sundrum (RS) type [4J 
can be achieved by considering as "bulk" the available space between the two parallel branes. 
Note, that in the case of the first RS model the size of the extra dimension is very short: in 
particular it is of the order of the Planck length. In such a case, gauge fields and fermions 
are not necessarily localized on the brane, see for example [9] and references therein. For 
our purposes in this article, we call this second scenario the Extended String /Brane- World 
scenario (ESBW). 

There are numerous generalizations of the above generic models, including, for instance, 
topological defects along the extra dimension(s), or higher-order curvature corrections, see 
for example Refs. [Tni[Ill[l2l[13l[Tail5lll6l[I71[T8l[T9l[20]. Also, there are models in 
which the standard model particles are localized on the brane dynamically, via a mechanism 
which is known as layer-phase mechanism, see for example Refs. [21], [22l [23l [2H [25l [26] and 
references there in. Moreover, an effective propagation of standard model particles in the 
bulk may characterize the so-called "fuzzy" or fluctuating-thick-brane- world scenarios [27] . 
according to which our brane world quantum fluctuates in the bulk. In such a case, there 
are uncertainties in the bulk position of the brane world, resulting in an "effectively" thick 
brane [28]. In such a situation, even the standard SSBW scenario would result in an "ef- 
fective" propagation of the standard model particles in the bulk. In such "effectively thick" 
brane models [28], the wavefunction of a standard model particle has an extent in the bulk 
region, the latter being defined as the "effective" thickness of the "fat brane" , of string (or 
Planck) size. 

In the previous part we presented a short review on brane world models, in what follows 
we focus to the so called asymmetrically warped brane models, and we give the essential 
references for the understanding of this paper. In particular, we will consider the following 
generic ansatz for the metric in five dimensions 

ds^ = -a\z)dt^ + (3^{z)dyi^ + -i^{z)dz^ (1) 

where z parameterizes the extra dimension, a{z) is the time warp factor and (3{z) is the 
three-dimensional-space warp factor |^. Models which are described by metrics of the form of 
Eq. ([1]), in which the time and three-dimensional-space warp factors are different, are often 
called asymmetrically warped brane models. In these models, although the induced metric 
on the brane (localized at ^ = 0) is Lorentz invariant upon considering the case a(0) = /3(0), 
the metric of Eq. ([T]) does not preserve 4D Lorentz invariance in the bulk since a{z) ^ (3{z) 
for 2; 7^ 0. 

Models with equal warp factors, such as the RS model [H [5] , have so far attracted the 
main attention, since 4D Lorentz invariance is assumed as a fundamental symmetry of nature. 
However, brane models with asymmetrically warped solutions, like that of Eq. ([T]) above, 
have also been constructed [221 EHl ED [321 [SS]- We note here that the static asymmetric 
solutions that are presented in Refs. [501 [SD 1321 [33] presuppose the existence of extra matter 
in the bulk [j, hence these models luck the simplicity of construction characterising the RS- 



^Although one can set 7 to one by a coordinate transformation in static models, this is not the case in 
more general cases with a time dependent 7(2, i), which are of interest in cosmology [29] . and which we make 
use here. It is in this sense we base our analysis here in the frame where the 5D metric has the form (jlj. 

^For example, the solution or Ref. |30j requires a constant electric field along the extra dimension, and a 
fine tuning similar to that of the RS model. 



model. 

Although symmetric |j brane models are the most important, as they incorporate exact 
4D Lorentz symmetry, at least in our opinion there is an argument in favor of asymmetrically 
warped space-times, presented in detail in Ref. [29]. Static solutions, like that of RS model, 
can be used approximately only for a short period of time around our cosmological epoch 
to- For larger periods the complete cosmological evolution must be considered. In such a 
cosmological context the metric depends on the comoving time coordinate t, and has the 
form: 

ds"^ = -a'^{z, t)dt^ + /^^(z, t)d^^ + 72(;z, t)dz^ (2) 

Note, that the above time depended solutions have in general different space and time warp 
factors. The corresponding static metric, which describes our universe in the short time 
period humans exist, which is very small compared to the cosmological time scale of the 
evolution of the Universe, is obtained if we set t = to'vn Eq. ([2]). According to this philosophy, 
in general, the metric of Eq. ([2]) evolves to an asymmetrically warped solution for t = to, or 
equivalently we expect that a{z,to) ^ /5(z, to)- Although, we have given arguments on the 
possibility that such vacua can indeed characterize cosmological brane models, nevertheless 
we cannot say with certainty that they cannot decay to standard Minkowski vacua. However, 
since they cannot be ruled out on theoretical grounds, at least presently, it worths imposing 
constraints on their parameters phenomenologically, and this is what we aim to do in this 
article. 

The question then arises as to how one can constrain or exclude/falsify brane models with 
asymmetric solutions of the form of Eq. ([T]), on account of present (or immediate-future) 
experimental bounds on (local) Lorentz symmetry violation in the bulk. In the SSBW case, 
where the bulk is completely unaccessible by the standard model particles, Lorentz violation 
signals can be observed only by bulk particles in the gravitational sector, like gravitons, or at 
most particles neutral under the standard model group, e.g. right-handed sterile neutrinos. 
Bulk fields can " see" the asymmetry between the warp factors in the extra dimension, whilst 
standard-model particles, which are rigidly "pinned" on the brane, can only "see" equal warp 
factors a(0) = (3{0). Gravity effects which could reveal 4D Lorentz violation are described 
in Ref. [31] (see also Refs. [Ml [35l [36]) where superluminous propagation of gravitons 
is possible for specific models with asym,m,etric solutions. However, since the detection of 
gravitons is still not an experimental fact, such Lorentz violations are still compatible with 
the current experiments, both terrestrial and astrophysical, probably awaiting the future 
detection of gravitational waves in order to be constrained significantly. 

However, in the ESBW scenario, where some or all of the standard model particles, are 
allowed to propagate in the bulk, such Lorentz-Invariance-violating effects can be bounded 
by high precision tests of Lorentz symmetry, since now, 4D Lorentz violation may be revealed 
even in the standard model sector. In this way, stringent restrictions to asymmetric models 
can be imposed by astrophysical observations and other high-energy experimental tests. 

In this paper, we aim to study the propagation of bulk photons in an asymmetrically 
warped metric background by solving the (classical) equations of motion for a 5D massless 
U(l) Gauge field in the curved background of Eq. ([T]). We shall demonstrate that the 



■^With the terminology symmetric solutions we mean here metrics with equal space and time warp factors, 
or equivalently metrics which possess S0(3,l) symmetry, like that of RS model. 



standard Lorentz invariant dispersion relation for 4D photons possesses nonlinear corrections, 
which lead to an energy-dependent speed of light on the brane. Specifically, we shall obtain 
a sub-luniinal refractive index for photons neff{uj) = 1 + cg uj"^, where u is the energy of 
the photon, and the factor cq is always positive and depends on the free parameters of 
the models. Finally, comparing these results with the data of the Magic experiment on a 
delayed arrival of highly energetic photons from the distant galaxy MkSOl [STl EB] we shall 
set concrete restrictions to our models. We shall also discuss briefly other astrophysical 
constraints on our models. 

2 5D U(l) Gauge fields in asymmetrically warped space- 
times 

In this section we study bulk photons in an asymmetrically warped metric background of 
the form ([1]). Although, our conclusions are valid for the quite general ansatz of Eq. ([1]), 
nevertheless it is convenient to demonstrate our results for specific classes of solutions to 
the 5D Einstein Equations. In what follows, we shall consider a slightly modified RS-metric 
which is asymmetrically warped, and it satisfies the bulk Einstein equations. We shall 
consider a class of solutions which are known as 5D AdS-Reissner-Nordstrom Black Hole 
Solutions [31]. These solutions are presented in detail in the following subsection, 2.1. In 
subsection 2.2 we show that is possible to write the 5D AdS-Reissner-Nordstrom Solution 
as a linearized perturbation around the RS metric. Then, it is not difficult to mimic the 
S^ /Z2 structure of the first RS-model |1]. However, the novel feature of our model is that the 
background metric deviates slightly from the standard RS-metric, and as a result there is a 
small violation of 4D lorentz symmetry in the bulk. In section 2.3 we study the propagation 
of bulk photons in the above-mentioned metric background. At this stage, we consider it as 
instructive to remind the reader that the extra dimension in the first RS-model is assumed 
to be compact and very small, of the order of the Planck length {ip ~ 10~^^cm). Hence, our 
assumption of bulk Gauge field is not in confiict with the current observations in particle 
accelerators, as we have also emphasized in the introduction. 

2.1 5D AdS-Reissner-Nordstrom Solutions 

We consider an action which includes 5D gravity, a negative cosmological constant A, plus 
a bulk U(l) Gauge field 



S = j (fx^ (^_i_(i?(5) _ 2A) - -^B"'''Bmn) + / d^X^J^j;;—^Cmatter, (3) 

where G5 is the five dimensional Newton constant, and Bmn = duHN — OnHm is the field 
strength of the U(l) Gauge field Hm-, with M, A^ = 0, 1 ... 5. Note, that this additional 
bulk Gauge field does not interact with charged matter on the brane, so it must not be 
confused with the usual electromagnetic field Am-, representing a bulk photon, which will be 
introduced in the next subsection. The four dimensional term in the action corresponds to 
matter fields localized on the brane, which is assumed located at r = tq, and are described 
by a perfect fluid with energy density p and pressure p. This brane term is necessary for 
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the solution of Eq. dH]) to satisfy the junction conditions on the brane, as we will explain in 
detail in the next section. The corresponding Einstein equations can be written aqj 



I Q(brane) I 

Gmn + AgAm = SvrGs | (5(r - ro) ^ '^ ' T^r"^'-),^^,^]^ + T^ 1 (4) 

/l^l 

where the energy momentum tensor for the U(l) Gauge field is: 

^MN = BmpBj^ — -qmnSpsB (5) 

For the metric of the black hole solution we make the ansatz 

ds^ = -h{r)de + rV^fiS^ + h{r)-^dr^ (6) 

where dT? = da'^+a'^dfl'^ is the metric of the spatial 3-sections, which in our case are assumed 
to have zero curvature, in agreement with the current astrophysical phenomenology, pointing 
towards spatial flatness of the observable Universe. Moreover, i is the AdS radius which is 
equal to y — f- -^y solving Einstein equations (jl]) we obtain: 

"(0 = ^-^ + ^ (7) 

where n is the mass (in units of the five dimensional Planck scale) and Q the charge of the 
5D AdS-Reissner-Nordstrom black hole. This, of course, presupposes the existence of extra 
bulk matter, namely a point-like source with mass /i and charge Q. Note that, in the case 
of nonzero charge Q, a non- vanishing component Bq,. of the bulk field-strength tensor Bmn'- 

is necessary in order for the solution to satisfy the pertinent Einstein-Maxwell equations. 

In the limit Q = 0, which corresponds to the 5D AdS-Schwarzschild Black Hole, the 
bulk vector field Hm becomes redundant. In this article we shall also examine this case, in 
connection with effects on photon propagation in 4D and, as we shall see, the presence of 
the extra charge in the Reissner-Nordstrom case does not affect qualitatively the conclusions 
regarding the order of magnitude of the asymmetry between the warp factors in order to get 
consistency with the current phenomenology. 

2.2 Construction of two-brane models and Null Energy Condition 

We will describe now how one can use the 5D AdS-Reissner-Nordstrom vacuum in order to 
construct brane models. As a first step, we place a brane at the position r = tq. We next 
assume that for r > Tq the 5D metric is given by Eq. (^, while for r < Tq the metric is given 
by Eq. (jSj) upon the replacement r <-^ Tq/t. The metric obtained this way is Z'2-symmetric 



*In the case of a second brane at the position t'q we need an additional 6{r — rg) term in equation Q. 



upon the replacement r ^-> Tq/t. The next step is to glue the two independent slices of the 
metric by including a perfect fluid energy momentum tensor of the form 

T^'' = Diag(-p,p,p,p) (9) 

and satisfying the Israel junction conditions on the brane, where p is the energy density and 
p is the pressure. The equation of state is parametrised as usual by: u = p/p. 

We remind the reader, that 5D black holes suffer from a singularity at r = 0, which may 
be naked or not, depending on the magnitudes of the parameters Q,p,i. In particular, if 
Q^ < ^p^i"^ the 5D black hole has two horizons [31j. In this case, a single brane construction 
is possible as a horizon can isolate the brane from the singularity. Contrary to the Randall- 
Sundrum (RS)-model, this class of brane models has the advantage that one does not need 
a fine tuning between the energy density p on the brane and the 5D cosmological constant. 
This self-tuning property of the AdS-Reissner-Nordstrom black holes serves as a possible 
solution to the cosmological constant problem [31] . 

However, if we demand the existence of a horizon, the matter we have to put on the 
brane violates the null energy condition p + p > (or tu > —1), see Ref [31]. This is a quite 
embarrassing situation, as we have to introduce unconventional sources like ghosts. In Refs. 
[Ml [35] the authors present a complete investigation on this issue including cases with: a) 
the presence of external fields in the bulk, b) the existence of higher-derivative terms in the 
action of Gauss-Bonnet type and c) models with the Z2 symmetry relaxed. The result of 
such analyses is a no-go theorem, according to which it is impossible to shield the singularity 
from the brane by a horizon, unless the null energy condition p + p > is violated. 

Nevertheless, in the current work we are not concerned with a solution to the cosmological 
constant problem. Our purpose is simply to construct brane models which incorporate a 
Lorentz symmetry violation due to the difference of the space and time warp factors. In this 
case we can isolate the singularity from the brane by putting a second brane at a position 
r = Tq (between the first brane and the singularity). Hence, we have a two- brane construction 
similar to the first RS-model, in which the two branes are located at the fixed points of an 
orbifold S1/Z2. Note, that because the Israel junction conditions must be satisfied also on the 
second brane, we have to assume an additional energy momentum tensor of a perfect fiuid, 
with energy density p', pressure p' and uj' = p'/p', localized on the second brane. We also 
consider the physically interesting case in which e = TQ/ro ~ 10~^ where the large hierarchy 
between the Planck and electroweak scale is satisfied. However, in the case of two-brane 
models fine-tuning relations between 4D energy densities and 5D dimensional cosmological 
constant are unavoidable. 

The question that arises at this point is wether two-brane constructions that respect the 
null energy condition on the branes do exist. An analysis towards this direction can be found 
in Ref. [36]. We have also performed independently a similar analysis, which confirms these 
results, but it is quite lengthy to be presented in this paper, as this would take one out of the 
main point of this work. For our purposes here, it suffices to state that construction of several 
classes of two-brane models in the 5D AdS-Reissner- Nordstrom vacuum is possible, without 
violation of the null energy condition, and hence two-brane constructions are acceptable. 



2.3 5D AdS-Reissner-Nordstrom Solution as a linearized pertur- 
bation around the Randall- Sundrum metric 

In order to achieve our aim, which is to write the 5D AdS-Reissner-Nordstrom solution as a 
linearized perturbation around the RS metric, we perform the following change of variables 
r — »• z{r) in Eq. IQ: 

r = roe"'' "^ , for z > 

r = roe''^ , for z < , (10) 

If, in addition, we make the rescaling x^ — > yx^ (yU = 0, . . . 3), we obtain: 

ds"^ = -a^{z)h{z)dt^ + a^{z)d^^ + h{z)-^dz^ (11) 

where a{z) = e"*^'^', and k = i~^ is the inverse AdS^ radius. For the function h{z) we obtain: 

h{z) = 1 - 6h{z), 6h{z) = /^e^'^l^l - ^e''\'\ (12) 

As we describe in detail in section 2.2, it is not difficult to construct two brane models. 
Now, the Z2 symmetry r <-;> r^/r if it expressed in the frame of the new parameter z, reads 
z — >■ —z. In addition, the positions of the branes which are located at tq and Tq = r^e in the 
original coordinate system, in the new coordinate system are determined by the equations 
z = and z = Trr^ correspondingly, where e = e'^'^^" and r^ is radius of the compact extra 
dimension. 

Finally, we assume that 1 5/^(2;) | -C 1 in the interval < z < nrc, or equivalently we 
adopt that 6h{z) is only a small perturbation around the RS-metric. This implies that ro, 
which is the radius that determines the position of the brane in the bulk, is (comparatively) 
a very large quantity. In particular, we have to satisfy both the following two inequalities 
rge^ > ^/J^(^ and r^e^ > Qi. 

2.4 Equation of motion for the Bulk Photon and Non-Standard 
Vacuum Refractive Index 

In this subsection we will study the case of a 5D massless f/(l) gauge boson A^ in the 
background of an asymmetrically warped solution of the form of Eq. ( ITT1) . We stress that 
the gauge field Aj^ must not be confused with the gauge field if at, introduced in the previous 
section. As we will see later, we will identify the zero mode of A^ with the standard four 
dimensional photon. On the other hand H^ is an additional bulk field which does not 
interact with the charged particles on the brane. The corresponding equation of motion for 
An reads: 

-^Om (V^g'^'^g^'F^s) = , (13) 

In the background metric of Eq. (iTTl) . 
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where we have considered the Coulomb gauge condition: 

V-i=0, Ao = 0, A, = 0, (15) 

which is suitable for the case of a Lorentz violating background. On setting in Eq. fll3p : 

Aj{x, z) = e'P-'^Xjiz), Pf. = (-W, p) (16) 
and keeping only the linear terms in 6h{z), we obtain 

- d, {a\z) [1 - 6hiz)] d^x} + {p' - [1 + Sh{z)] uj'}x = (17) 

where for brevity we have dropped the index j from x- Note that the spectrum of Eq. (TT71) 
is discrete, due to the orbifold boundary conditions [1], x'(0) = and x'(^^c) = (where 
the prime denotes a ^-derivative). Hence we add from now on an index n (n = 0, 1, 2, 3, ..), 
denoting the nth eigenstate of the equation. 

Upon applying the formalism of time-independent perturbation theory, Eq. flT7|) can be 
written as 

(Ho + AH)xn = ml Xn, m^ = to^ - p"" . (18) 

The hermitian operator 

Ho = -d, {a'iz)d,) , (19) 

which is of zeroth order in 6h{z), corresponds to the unperturbed Hamiltonian, while the 
hermitian operator 

AH = -6h{z) u^ + d,{a^{z)5h{z)d,) (20) 

which is linear in 5h{z), corresponds to the perturbation Hamiltonian. 

The quantity AH is indeed a small perturbation when compared to Hq, since the relative 
discrepancy 6h{z) between the space and time warp factors is assumed to be very small 
(|5/i(z)| << 1) in the interval [0,7rrc]. 

We need to find a complete set of eigenf unctions Xn^ (^ = 0, 1, 2, ...) of the unperturbed 
Schrodinger equation 

^0X1°^ = K0'^"^ (21) 

It can be seen straightforwardly that for n = 0, 

^r = , (22) 

is an eigenvalue of Eq. fl2Tl) . or equivalently Eq. fl2Tl) possesses a zero mode, which corre- 
sponds to the 4D photon. The corresponding wave function is constant in the interval [0, vrrj. 
For our purposes here we shall omit the detailed derivation of the eigenvalue spectrum of 
(!2T|) . as this has been done previously in Refs. [391 HO], where we refer the interested reader 
for details. It suffices to mention that the nonzero eigenvalues are: 

(^')=Xnke-''-'%n = l,2,3,... (23) 



m, 



where Xn are the roots of the zeroth order Bessel function Jo(a;„) = 0. On adopting kvc ~ 12, 
which is the standard choice in order to connect electroweak (ew) and Planck scales {Mp = 



-^k-KTc 



nrip 



in a RS framework [1], one obtains that: 



m. 



f~TeV, n = l,2,.... (24) 



The corresponding eigenfunctions are: 



xr = ^, iVo = -L (25) 

No y/W2 

Xf = Tre'Vi(^e^^), iV„ = ^Ji(x„) (26) 



Nn '' k " " V2k 

where the coefficients A'^o, iV„ are defined by the normahzation condition: 

X^:Hz)x^^\z)dz = 6^n (27) 



irrc 
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The lowest energy eigenvalue ml of Eq. (TTS!) can be expanded as: 

>2 / n^\2 / /9^^2 



^^ = Kt + Kt + Kt + - (2j 



where the zeroth order correction (zero mode) is ?Tig = (c.f. (!22l) . The first and second 
order corrections are obtained by the formulae of time-independent perturbation theory: 

(4^^ = (xi°V^lxr) (29) 

. (2).. ^ y. KxriAg|xi°))P .30. 

Upon using Eqs. ([28]), ([29]) and ([30]), we obtain: 

tu^ — p^ = — ac uj'^ — bo u;^ (31) 
with 



og 



rfz(xr(^)) 5M^) (32) 



^G 



E ^^^ (prf^xr(^)xi°n^)'5M^)) (33) 



Note, that the derivative term in (!20i) does not contribute to a^ and be, as it acts on a 
constant eigenfunction (zero mode). The parameters ac and bo are independent from the 
(photon) energy u, and 6^ is always positive. Eq. (|3T]) then gives: 



|p|=a;^(l + aG) + 6G^2 (34) 

According to Eq. (IMl) . the phase velocity of light is equal to 

UJ 1 

Vph = ri = I (35) 

IpI J{l + aG)+bGu;^ 



or equivalently, for bc^^ << 1: 



Vph - I — :t^ (36) 

'l + aa) 2(l + aG)5 



where uj is the energy of the photon. 

On the other hand, for the photon's group velocity we have: 



_d^_ v/(l + aG)+feG^^ 



and the corresponding formula for 6g << 1 is 



Var — = T^ 38 

'[l + ac) 2(1 + aG)t 



One can define the constant velocity of light in standard vacuo as the low energy limit 
{uj -^ 0) of the phase velocity of Eq. (l36l) 

1 , , 

Clight = I (39) 

[l + ac) 



This definition is reasonable, if we take into account that the velocity of light is measured in 
practice only for low energy photons. We need very high energies, in our case of order 10^ 
TeV or higher, in order for the contribution of the energy dependent term in Eq. ( I36l) to be 
significant, exceeding the accuracy with which the velocity of light is presently known. 

From Eqs. fl36|) and fl39l) we then obtain an effective suhluminal refractive index for the 
non-standard vacuum provided by our brane-world constructions: 

n (^) = ^^1^ = 1 + ^^ y . (40) 

"^ ' Vph 2(1 + aG) ^ ' 

This is the main result of this paper. An issue, we would like to emphasize, is that the phase 
and group velocities (Eqs. (136|1 and (138!) respectively), as well as the effective refractive index 
of Eq. (HOj) . depend quadratically on the photon energy oj. Moreover, we note that equation 
(HUj) is a perturbative result which is vahd only for energies uj « 1/v^6g- For larger energies 
(for example in the case of ultra high energies cosmic rays uj ~ 10^° eV), Eq. (HOl) is not 
valid and a full nonperturbative treatment is necessary. In this limit the subluminal nature 
of the refractive index may be lost. 

According to the refractive index (HOj) . photons with different energies propagate with 
different velocities. Due to the sub-luminal nature of the refractive index there will be 
delays in the arrival times of more energetic photons, as compared with the lower-energy 
ones, provided of course that simultaneous emission of the various modes is assumed at the 
source. Hence, we will observe a time lag of the arrival times of photons which were emitted 
simultaneously by a remote astrophysical source, with the lower energy photons arriving first. 
It is worth noticing that the effective refractive index (HOl) is independent on the polarization 
of photons, or equivalently our model does not predict birefringence. If we had included 
higher order corrections to the electromagnetic lagrangian (e.g. in the form of a Born Infeld 
lagrangian) birefringence may occur. In the next section (section 3) we compare our result 
with the data of Magic experiment in order to set bounds to our model. 
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3 Setting bounds: Comparison with MAGIC and other 
Data 

In the previous section we examined the case of bulk photons in asymmetrically warped brane 
models and we derived an equation (Eq. (HOj) ) which implies a nontrivial refractive index 
for 4D photons. In this section we aim to use this modified quadratic dispersion relation 
in order to set bounds to our model, via a direct comparison with the experimental data of 
Magic. We have used the MAGIC Telescope observation of photons from MkSOl [371 EH], 
because this is the only case at present of an observed effect on a delayed photon arrival, 
with the higher energy photons arriving later. The constraints we derive, in this section, 
for the parameters of our model, are the most stringent on quadratically modified photon 
dispersion relations, as far as direct photon observations are concerned. Note, that there are 
other experiments which set bounds on quadratically modified photon dispersion relations, 
however they do not give more stringent constraints than Magic effect. For example in 
subsection 3.2, we compare with the recent data from H.E.S.S. Telescope [H], and as we see 
the magic constraints are the most stringent. Finally, in subsection 3.2, we discuss briefly 
possible constraints to our model which are not based on the quadratic dispersion relation 
of Eq. ( I40l) . including the cases of the fermions and the ultra high energy cosmic rays. 



3.1 Comparison with MAGIC Data 

We wish now to compare the time delays of the more energetic photons, implied by the 
dispersion relation (HOj) . with experimentally available data. In view of the quadratic sup- 
pression of the non trivial vacuum effects with the Planck scale, it becomes evident that the 
higher the energy range of the photons detected, the higher the sensitivity of the pertinent 
experiment to such effects. Recently, experimental data from observations of the MAGIC 
Telescope [3^ on photon energies in the TeV range have become available. It therefore 
makes sense to compare the time delays predicted in our models, due to (HOj) . with such 
data, thereby imposing concrete restrictions on the free parameters of the models. 

For completeness let us first review briefly the relevant observations [371 [38] . MAGIC is 
an imaging atmospheric Cherenkov telescope which can detect very high energy (0.1 TeV-30 
TeV) electromagnetic particles, in particular gamma rays. Photons with very high energy 
(VHE) are produced from conversion of gravitational energy at astrophysical distances from 
Earth, when, say, a very massive rotating star is collapsing to form a supermassive black 
hole. Astrophysical objects like this are called blazars and are active galactic nuclei (AGN). 

The observations of MAGIC during a flare (which lasted twenty minutes) of the rela- 
tively nearby (red-shift z ~ 0.03) blazar of Markarian (Mk) 501 on July 9 (2005), indicated 
a 4 ± 1 min time delay between the peaks of the time profile envelops of photons with 
energies smaller than 0.25 TeV and photons with energies larger than 1.2 TeV. Possible 
interpretations of such delays of the more energetic photons have already been proposed. 
Conventional (astro)physics at the source may be responsible for the delayed emission of 
the more energetic photons, as a result, for instance, of some non-trivial versions of the 
Synchrotron Self Compton (SSC) mechanism [37]. It should be noted at this stage that the 
standard SSC mechanism, usually believed responsible for the production of VHE photons 
in other AGN, such as Crab Nebula, fails [37] to explain the results of MAGIC, as a relative 
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inefficient acceleration is needed in order to explain the (9(min) time delay. Modified SSC 
models have been proposed in this respect [12], but the situation is not conclusive. This 
prompted speculations that new fundamental physics, most likely quantum-gravity effects, 
during propagation of photons from the source till observation, may be responsible for in- 
ducing the observed delays, as a result of an effective refractive index for the vacuum, see 
for example Ref. |43j and references therein. 

Exploring the fact that MAGIC had the ability of observing individual photons, a nu- 
merical analysis on the relevant experimental data has been performed in [38] , which aimed 
at the reconstruction of the original electromagnetic pulse by maximizing its energy upon 
the assumption of a sub-luminal vacuum refractive index with either linear or quadratic 
quantum-gravity-scale suppression: 

n,ff{uj) = l+(-^] , n = l,2 (41) 

The analysis in ^38j resulted in the following values for the quantum-gravity mass scale at 
95 % C.L. 

Mqg(i) - 0.21 10^^ GeV, Mqg(2) ^ 0.26 10^^ GeV (42) 

We must emphasize at this point that many Quantum Gravity Models seem to predict modi- 
fied dispersion relations for probes induced by vacuum refractive index effects, which appear 
to be different for each quantum gravity approach, not only as far as the order of suppres- 
sion by the quantum gravity scale is concerned, but also its super- or sub-luminal nature. 
Some models, for instance, entail birefringence effects, which can be severely constrained by 
astrophysical measurements [l^. There are also alternative models, see for example Refs. 

Ha IS]. 

In |17], a non-perturbative mechanism for the observed time delays has been proposed, 
based on stringy uncertainty principles within the framework of a string/brane theory model 
of space-time foam. The model entails a brane world crossing regions in bulk space time 
punctured by point-like D-brane defects (D-particles). As the brane world moves in the bulk, 
populations of D-particles cross the brane, interact with photons, which are attached on the 
brane world (SSBW scenario), and thus affect their propagation. The important feature of 
the stringy uncertainty delay mechanism is that the induced delays are proportional to a 
single power of the photon energy u, thus being linearly suppressed by the string mass scale, 
playing the quantum gravity scale in this model. It is important to notice that the mechanism 
of [17] does not entail any modification of the local dispersion relations of photons. In fact 
the induced delays are also independent of photon polarization, so there are no gravitational 
birefringence effects. It is also important to notice that, in view of the electric charge 
conservation, the D-particle-foam defects can interact non trivially only with photons or at 
most electrically neutral particles and no charged ones, such as electrons, to which the foam 
looks transparent. 

As we have showed in the previous section, a sub-luminal vacuum refractive index may 
characterize asymmetrically warped brane-world models, assuming propagation of photons 
in the bulk. Our model, however, predicts a refractive index with quadratic dependence on 
energy, c.f. fj40|) . Comparing Eqs. fHOj) with PTI) . we obtain: 

^"^ < ^QG(2) (43) 



2(1 + aG 
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The parameters a^ and bo can be computed numerically by means of Eqs. ([32D and 
( !33|) . Note, that the deviation 5h{z), the eigenvalues m!^^ and the eigenfunctions Xn^ cire 
known analytically, see Eqs. (TT2|) . fl23l) . fl25|) and (126!) above. We have computed numerically 
the parameter ba for two exemplary cases: a) for an AdS-Reissner-Nordstrom Solution (TT2|) . 
and b) for an AdS-Schwarzschild Solution (obtained by setting Q = in Eq. (1121) ) . In the 
former case we assume for simplicity that /i and Q^/tq are of the same order of magnitude 
(for details see Ref. [36] or the discussion in section 2.2, see also Eq. fl53|) below), and hence 
we can ignore the contribution of the e^'^^ term in Eq. flT2l) . In particular we obtain: 

bo ^ 2.95 ((5/i)^ TeV"^, Ads - Reissner - Nordstrom (44) 

be ^ 10 {Sh)^ TeV-^ Ads - Schwarzschild , (45) 

where {Sh) is defined as the average value: 

{Sh) = / Sh{z)dz . (46) 



vrrc JO 

We will use {6h) in order to estimate the degree of violation of Lorentz symmetry in our 
models. Taking into account Eqs. P2l) . fB3|) . fjl^ and fH5l) we find the constraints: 

I {6h) |< 1.4 10"^ Ads - Reissner - Nordstrom (47) 

I {6h) \< 0.75 10"^ Ads - Schwarzschild . (48) 

The small values we obtain are consistent with the weak nature of (Sh), as required by 
treating it as a perturbation. Also we observe that the values for the average deviation 
of Eqs. (jHj) and (HTIl are of the same order of magnitude for both Ads- Schwarzschild and 
Ads-Reissner-Nordstrom solutions. In the above analysis we ignored the effects due to the 
Universe expansion, since the latter do not affect the order of magnitude of the above bounds 
due to the small red shift {z ~ 0.03) of Mk501 we restrict our discussion in this section. The 
inclusion of such effects, which are essential for larger redshifts, is straightforward [38j and 
does not present any conceptual difficulty. 

Note, that although the parameter a^ is not important for our analysis, it is crucial when 
we have to make comparisons with the velocities of other particles. If we take into account 
Eq. (1521) we see that ac = {Sh). Hence, the parameter ac is constrained via the equation 

I flG |< 10"^ (49) 

This summarizes the constraints to the asymmetrically-warped brane models with bulk pho- 
tons using the data of the MAGIC experiment. 

3.2 Comparison with other data 

However it is worth to be mentioned that another Telescope, H.E.S.S., looking at a flare 
from another celestial object, the Active Galaxy PKS 2155-304 at a much larger distance 
(red-shift z = 0.116), and emitting photons in the few-TeV energy range with a much higher 
(unprecedented) statistics, has not claimed an observed delayed arrival of the most energetic 
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photons. The associated hmit in the quadratic refractive index from such observations 
was (41j : 

Id < 7.3 X 10^^ , n^ffiuj) =(l + C^] (50) 



with Mp = 1.22 X 10^^ GeV, the conventional Planck mass scale in the notation of Ref. ^T] . 
Using fHOl) we can translate these bounds to bounds for the parameters of our model, in an 
entirely analogous way with the MAGIC case fl47|48l) studied above. Because of the larger 
distance from Earth, however, of the AGN PKS 2155-304 as compared to Mk501, the effects 
of the expansion of the Universe [121 EH] niust be taken into account in this case, when 
fitting the results to the observed arrival times of photons. We shall not do this in this work, 
since the bounds obtained from the H.E.S.S. measurements fl50l) are not better from the 
ones obtained from MAGIC, discussed above. Indeed, in the notation fl50|) the lower bound 
of the quantum gravity scale from these null measurements (as far as delayed arrivals are 
concerned), assuming quadratic modifications in the photon dispersion relations, is [H] 



Mqg(2) > Mp/ICI^/^ ~ 1.4 X 10^ GeV , (51) 

which is worse than the corresponding MAGIC bound ( H2l) . and hence for our purposes 
here the MAGIC bound is the most stringent one based on direct high-energy astrophysical 
photon time-of-flight measurements. 

However, in the models examined in this paper the modified dispersion relations are 
expected to characterize also other massive particles that may be allowed to propagate in 
the bulk. One may in this way impose more severe restrictions, for instance, by comparing 
the velocity of bulk photons with that of other massive particles, such as massive fermions. 
For such a comparison, additional considerations are necessary: a) the fermions are localized 
in one of the two branes or b) the fermions can propagate in the bulk, a case which has 
not been analyzed in this paper. If, for example, one finds that fermions move faster than 
light, then severe additional restrictions have to be imposed, due to the emission of vacuum 
Cherenkov radiation, for details see Ref. |11] . However an investigation toward this direction 
is a topic of a completely new work. 

Moreover, in case a photon exhibits modified dispersion relations, as advocated here, then 
scattering of Highly Energetic photons (with energies higher than 10^^ eV) of soft photons 
of the Cosmic Microwave Background (7he + 7cmb -^ e~^e~) may lead to entirely different 
physics than in the Lorentz-invariant case, as this affects the energy thresholds for pair 
production. This may lead to severe constraints on the relevant models [IB]. Naively, one 
could use the bounds, which are presented in Ref. [H] , to set limits in the parameters of our 
model. However, the analysis of such very high energy photons requires a nonperturbative 
analysis going beyond the approximations made in the present paper, and is left for future 
investigation. 

4 Comparison with the velocity of gravitons 

Finally, before closing we would like to compare our results for bulk photons with the velocity 
of gravitons, which are also allowed to propagate in the bulk. In Refs. pT| [M] the following 
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relation is found: 

12/72 



(bulk) _ , V_f_^^-4 _ 0^-2^ 

graviton — -L "T ^^g 



^graviton ~ ^ ~^ a & \^ ) v^ / 



where we have adopted the position of Ref . [36] , in the case where the null energy condition 
on the brane is not violated, ioi w = —1: 

/^ = oTT (53) 

On the other hand, as discussed here, bulk photons propagate with a velocity given by Eq. 
(139|) . or equivalently: 

«> = l-^. {ih)^aa (54) 

(«»> = 1^ i'^A^-* - 1) - '?^((-' - 1)) . « = e-"- ~ 10-"= (55) 



where 



2k7rrc \2rg 3r\ 

Taking into account Eq. (155]1 we find 






Then, the difference reads: 



^{bulk) (bulk) _, <5 ^ 



^Hlg) - flight araviton — ^2k7Tr r^ 

We observe that if photons propagate in the bulk then gravitons are not superluminous. 
However, we can not use this result to set constraints via gravitational Cherenkov radiation, 
since the velocity of bulk photons is in general different from the velocities of fermions, which 
we have not computed in this work. For more details on this topic the interested reader is 
referred to Refs. [36l HO] . 



5 Conclusions 

In this article we have examined non-standard Lorentz- violating vacua, induced in the brane- 
world framework due to asymmetrically warped bulk space times. We have considered the 
effects of such an asymmetric warping on four- dimensional photons, with wavefunctions that 
extend into the bulk. These assumptions lead to an effective refractive index of the 4D non- 
standard vacuum arising in this context, which scales quadratically with the quantum gravity 
(length) scale and photon energy. By comparing our results with data from the MAGIC 
experiment, pointing towards a delayed arrival time of more energetic photons from MkSOl, 
as compared to lower-energy ones, we have derived bounds on the size of the perturbations 
around the classical solutions, or order 10~^ (c.f. fl47ll48p ). This result seems rather generic, 
independent of the particular type of the assumed bulk solution. 

In order to set bounds, we have used the MAGIC Telescope observation of photons from 
MkSOl, because this is the only case at present of an observed effect on a delayed photon 
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arrival, with the higher energy photons arriving later. The constraints we derive, in section 
3, are the most stringent on quadratically modified photon dispersion relations, as far as 
direct photon observations are concerned. We would like to mention that there are also 
other experiments which set bounds on quadratically modified photon dispersion relations. 
In particular, in subsection 3.2, we compare with the recent data from H.E.S.S. Telescope 
[H], and we see that the magic constraints are the most stringent. Finally, in subsection 3.2, 
we discuss briefiy possible constraints to our model which are not based on the quadratic 
dispersion relation of Eq. fHOl) . including the cases of fermions and ultra high energy cosmic 
rays. However, an analysis towards these directions is quite extended to be included in this 
work, and is left for further investigation. 
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